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ABSTRACT 

Hypervelocity stars have been recently discovered in the outskirts of galaxies, such as the unbound 
star in the Milky Way halo, or the three anomalously fast intracluster planetary nebulae (ICPNe) in the 
Virgo Cluster. These may have been ejected by close 3-body interactions with a binary supermassive 
black hole (SMBBH), where a star which passes within the semimajor axis of the SMBBH can receive 
enough energy to eject it from the system. Stars ejected by SMBBHs may form a significant sub¬ 
population with very different kinematics and mean metallicity than the bulk of the intracluster stars. 

The number, kinematics, and orientation of the ejected stars may constrain the mass ratio, semimajor 
axis, and even the orbital plane of the SMBBH. We investigate the evolution of the ejected debris 
from a SMBBH within a clumpy and time-dependent cluster potential using a high resolution, self- 
consistent cosmological N-body simulation of a galaxy cluster. We show that the predicted number 
and kinematic signature of the fast Virgo ICPNe is consistent with 3-body scattering by a SMBBH 
with a mass ratio 10 : 1 at the center of M87. 

Subject headings: clusters, supermassive black holes, galaxies, n-body simulations 


1. INTRODUCTION 

A significant fraction of the stellar component of a 
galaxy cluster is not confined to any galaxy. This intra¬ 
cluster light (ICL) has been identified via planetary neb¬ 
ulae (PNe) (Feldmeier et al. 2004; Aguerri et al. 2005), 
Red Giant Branch (RGB) stars (Durrell et al. 2002), 
and ultra-deep surface photometry (e.g., Feldmeier et al. 
2002, 2004; Mihos et al. 2005). It is commonly thought 
that intracluster stars are stripped from galaxies as the 
cluster assembles (Merritt 1984), via high speed galaxy 
encounters (Moore et al 1996), interactions with the clus¬ 
ter potential (Byrd & Valtonen 1990), or by tidal strip¬ 
ping in infalling groups (Mihos 2004). These processes 
will generate a debris field that is highly inhomogeneous, 
with distinctly non-Gaussian velocities, reflecting an un¬ 
relaxed intracluster population (Napolitano et al. 2003). 

The continued search for intracluster PNe (IGPNe) 
has revealed several in the Virgo Gluster with radial ve¬ 
locities that are extremely rapid compared both to its 
nearest galaxy and to the background cluster. In one 
pointing ^ 65 kpc from M87, Arnaboldi et al. (2004) 
found that of 15 ICPNe radial velocities, only 12 were 
consistent with M87’s stellar velocity dispersion profile 
(Romanowsky & Kochanek 2001) and systemic velocity. 
Of the three remaining ICPNe, two had velocities which 
were offset ^ —900 km/sec from the systemic velocity of 
M87 and one had a ~ -1-1500 km/sec difference. 

While these ICPNe could be associated with an un¬ 
relaxed tidally-stripped intracluster stellar population, 
there is an alternative explanation. In this paper, we 
explore the possibility that the high velocity ICPNe 
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were ejected after interacting with a supermassive bi¬ 
nary black hole (SMBBH) at the center of M87. Super¬ 
massive black holes (SMBH) are a standard component 
of elliptical galaxies and spiral bulges (e.g., Richstone 
et al. 1998). When galaxies evolve through hierarchi¬ 
cal merging, each SMBH participates and at some point 
forms a hard binary. This implies that many galaxies 
host SMBBHs. Perhaps the best evidence for this is the 
double X-ray bright nucleii in NGC 6240 (Komossa et al. 
2003). 

During the hard binary phase, stars that pass between 
the two SMBHs siphon energy from the binary’s orbit via 
3-body scattering, and most are ejected with hyperveloc¬ 
ities (Hills 1988; Yu & Tremaine 2003). These SMBBH- 
driven ejecta may constitute a different ICL population 
than stars stripped by dynamical interactions between 
galaxies; they are more likely to be metal rich, and, as 
we will show, have distinct kinematics and spatial struc¬ 
ture. This makes the two populations easy to separate. 

We investigate the structure and kinematics of the de¬ 
bris ejected via 3-body scattering within a high resolu¬ 
tion N-body simulation of a cluster potential. We model 
the SMBBH ejection velocities and study how well the 
ejected population remains kinematically distinct within 
a clumpy and still evolving cluster potential. We com¬ 
pare the expected mass fractions of the ICL generated 
by 3-body interactions with more traditional tidal inter¬ 
actions. According to our models, it is plausible that the 
fast ICPNe were ejected after interacting with a SMBBH 
in M87. We discuss how this theory can be tested with 
spectroscopic observations. 

2. GENERATING HYPERVELOCITY STARS 

Galaxy mergers provide an impetus for SMBHs to meet 
and form a bound system. During a galaxy merger, each 
SMBH sinks to the center of the new galaxy potential due 
to dynamical friction, and eventually becomes bound as 
a SMBBH. Dynamical friction then continues to shrink 
the orbit until the binary is hard (i.e., the separation be- 
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tween each SMBH, obbh, is such that the system tends 
to lose energy during stellar encounters). Thereafter, fur¬ 
ther decay is mediated by 3-body scattering with the am¬ 
bient stellar background until the SMBBH becomes so 
close that the orbit can lose energy via gravitational ra¬ 
diation. It will thereafter presumably coalesce. The final 
stages of SMBBH coalescence emit so much gravitational 
radiation that they are extremely likely to be detected 
by the Laser Interferometer Space Antenna (LISA), a 
planned NASA mission to detect gravitational waves, set 
to launch in the next decade. 

When black holes form a hard binary, any star that 
passes nearby takes energy from the system and changes 
the star’s velocity (see Yu 2002): 
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where Mi and M 2 are the black hole masses, and ubbh = 
GM2/4CTg, where ctq is the central velocity dispersion. 
Since 3-body scattering preserves the z-component of the 
angular momentum, the ejected debris is confined to a 
torus with a major axis aligned with the instantaneous 
SMBBH orbital plane (Zier & Biermann 2001) . The 
kinetic energy exchange takes place over a timescale: 


6t ~ 


G(Mi -h M 2 ) 


3 

BBH 


( 2 ) 


To model this process in M87, we need to choose ap¬ 
propriate parameters for a SMBBH at the galaxy center. 
This is currently ill-constrained. If the jet in M87 were 
a long-lived structure, we could estimate the masses and 
separation of the binary by shape of jet, since a large 
secondary mass could induce periodic structure. How¬ 
ever, the jet’s length, along with a reasonable choice 
for the bulk gas velocity, implies that M87’s jet may 
only be 2-3 Myr old. Thus, it provides no real con¬ 
straint on the component masses. We therefore set the 
total mass of the SMBBH to 3.3 x 10®M© (Harms et al. 
1994; Macchetto et al. 1997 ). Semi-analytic models of 
the assembly and growth of SMBHs using an extended 
Press-Schechter formalism indicate that the most proba¬ 
ble mass ratio for coalescing SMBBHs at low redshift is 
M 1 /M 2 = 10 (e.g., Sesana et al. 2004). Hence, we set the 
ratio of the two SMBHs to be Mi/M 2 = 10. If we assume 
the binary is hard, then the separation is obbh =2.1 
pc, and Sv ~ 1350 km/sec. A full loss cone then im¬ 
plies a SMBBH coalescence timescale of O(IO^) years, in 
which case the SMBHs may already have merged (and 
may even be associated with the jet). However, less op¬ 
timistic assnmptions about the loss cone reservoir yield 
a much longer coalescence timescale, suggesting that the 
SMBBH may be observable. 


3. SIMULATING EJECTA KINEMATICS 

To compare the kinematics of the SMBBH-ejected 
stars to the stellar kinematics within a complex cluster 
potential, we use a high-resolution N-body galaxy cluster 
simulation as a testbed for our analytic 3-body ejection 
scheme (Mihos et al. in prep). To create this simu¬ 
lation, a 50^ Mpc^ A=0.7, Hm= 0.3 cosmological dark 
matter simulation was run from 2=50 to 2=0, at which 


point a collapsed cluster with mass ~ 10^'^ M© was cho¬ 
sen to re-simulate at higher resolution. At 2 = 2 , indi¬ 
vidual halos destined to end up within the 2=0 cluster 
were identified; halos more massive than 10% Milky Way 
mass were excised from the simulation and high resolu¬ 
tion collisionless galaxy models inserted in their place. 
To do this, we used the “halo occupancy distribution” 
formalism (Berlind & Weinberg 2002) wherein massive 
halos are more likely to have multiple galaxies inserted. 
Both spiral and elliptical galaxy models were used. Spi¬ 
rals consisted of an exponential disk plus a Hernquist 
(1990) bulge (with bulge:disk ratio of 1:5), while ellipti¬ 
cals were a pure Hernquist (1990) model. Both galaxy 
types were embedded in high resolution isothermal dark 
halos which blended smoothly into the background dark 
matter distribution. The galaxy models were scaled in 
mass to match the halo mass they replaced, scaled in 
size by M^/®, and scaled in velocity by M^/^. In total, 
I2I high resolution galaxy models were inserted into 80 
dark halos; the simulation consisted of approximately 10 
million particles, 5.4 million of which represented stars 
with a gravitational smoothing length of 0.28 kpc. Ini¬ 
tialized at 2 = 2 , the simulation was evolved to 2=0 using 
the N-body code GADGET (Springel et al. 2001). Fig¬ 
ure la shows the simulation at 2=0 assuming a stellar 
mass-to-light ratio of 5. 

To compare the characteristics of tidally-released stars 
with stars ejected via 3-body interactions, we must first 
identify the tidal ICL. Definitions of IGL abound, based 
on binding energy, morphology, or surface brightness. 
For our illustrative purposes, we simply identify particles 
as ICL if their local density falls below a threshold, set so 
that when the cluster is initialized at 2=2, the ICL frac¬ 
tion is essentially zero, and at 2=0 the ICL particles no 
longer trace the galaxy morphologies. As a consistency 
check, using this criterion we arrive at an ICL fraction of 
13%, similar to that inferred for Virgo from observations 
of ICPNe (Feldmeier et al. 2004). 

Even with such a high resolution N-body simulation, 
no particle would pass within obbh of our M87-analogue 
nucleus. We therefore developed a Monte Carlo approach 
to mimic the interaction between stars and a SMBBH. 
We pinned a spherical region with radius obbh to the 
center of the most massive cluster galaxy (with a mass 
comparable to M87), and embedded the cluster with 8000 
particles. We then gave these particles a small isotropic 
velocity dispersion so that they were destined to pass 
within Obbh rather quickly. Each particle received the 
velocity kick of equation 1 whenever it passed within 
Obbh of the galaxy center. The direction of this kick was 
planar, to ensure the debris would form a torus. Oth¬ 
erwise, all particles moved within the substructure-rich 
potential. This simulation was followed for 100 Myr. 

Figure lb and Ic depict the inner 300 kpc of the most 
massive galaxy. Here the torus of ejected PNe is plainly 
visible and the radial velocities are typically much faster 
than the ambient ICL. 


4. TOTAL MASS EJECTED BY SMBBHS 


We can estimate the total fraction of ICL ejected from 
M87 by a SMBBH during the hard binary phase: 


Mej = J(Mi -h M2)ln 
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Fig. 1. — Left: Surface brightness distribution of the inner 5 Mpc of the simulation. Colors represent the surface brightness in the 
V band. Middle: The projected position versus line-of-sight velocity for the inner 300 kpc of the simulation. The blue points are the 
particles ejected by the mock SMBBH at the center of the most massive galaxy. The ejected debris is clearly distinguishable from both 
the tidally-stripped ICL (green) and the galaxy itself (red). Right: Projection of the central 300 kpc around the largest galaxy with the 
observer tiled 50° above the orbital plane of the SMBBH. The torus thickness is related to the timescale over which the debris is ejected. 


where J is a dimensionless mass ejection rate. For equal- 
mass circular binaries, J ~ 0.5 (Milosavljevic & Merritt 
2001). Quinlan (1996) argues that J drops as M 1 /M 2 in¬ 
creases, and varies with the orbital speed of the SMBBH. 
We ignore the variation of J during the binary’s evolu¬ 
tion, setting J = 0.3. The separation where gravitational 
radiation dominates, Ugw is: 

_ ( 256 

agw - ^ 

where p is the central stellar density, c is the speed of light 
and H is the hardening rate, set by restricted 3-body 
experiments and direct N-body simulations. We adopt 
H — IQ \Ti the estimate that follows (Quinlan 1996), and 
use Faber et al. (1997) to set M87’s p and do- Given 
these constraints, we find that ~ 6 X 10® M 0 should be 
ejected by a 10:1 SMBBH in M87 during the hard phase. 
Although this is a significant amount of mass, it is only 
about 1.5% of M87’s total mass within 10 kpc (Nulsen 
& Bohringer 1995). Assuming M/L = 5, and that M87’s 
V-band luminosity is 7.7 x lO^^L©, nearly 2% of M87’s 
light could ejected during one 10:1 SMBBH merger. 

Scaling this result from M87 to the cluster as a whole 
requires a model for the SMBBH merger history of each 
Virgo galaxy, along with SMBH masses and mass ra¬ 
tios as the cluster assembles. To obtain a rough esti¬ 
mate, however, we can assume that each Virgo cluster 
galaxy hosted one 10:1 SMBBH merger in its lifetime, 
and ejected debris during the merger with an efficiency 
that depends on the black hole mass, M,. In reality, only 
a fraction of Virgo cluster galaxies would host SMBBHs 
mergers, yet each of these would likely have undergone 
several mergers hy z = 0. Our simple model averages 
over the galaxy ensemble and avoids the problem of iden¬ 
tifying mergers over a Hubble time. This reduces the 
problem to simply counting the appropriate galaxies and 
assigning a SMBBH mass to each one. 

Using the observed Virgo luminosity function 
(Sandage, Binggeli, & Tammann 1985), we take a census 
of the number of galaxies per absolute blue magnitude, 
between —22 < Mb < —12, as well as the variation 


of morphological type in each magnitude bin. Since 
SMBHs are thought to exist in every elliptical and spiral 
galaxy bulge, we include only spirals and ellipticals in 
our analysis. We model all spirals after the Milky Way, 
with the same mass, ao, and bulge to disk ratio of 1:3; 
our elliptical galaxies are modeled using Fundamental 
Plane parameters (Dressier et al. 1987). Naturally, 
spiral galaxies are more varied than the Milky Way, but 
as we will see, the choice of our spiral model matters 
very little for this exercise. 

Now that we have the number of SMBBH-embedded 
galaxies, as well as o-q for each, we determine the SMBBH 
mass in each galaxy using the M, — a relation (Gebhardt 
et al 2000; Ferrarese & Merritt 2000). Unfortunately, 
there is some disagreement on the precise slope of the 
relation, so we adopt M, oc (Tq (e.g., Tremaine et al. 
2004). Since we have the SMBBH mass and o-Qj we also 
have ubsh and Ogw, and hence Mej for each galaxy. 

Folding in these various relationships and summing 
over the luminosity function, we obtain a total SMBBH- 
ejecta luminosity ~ 2% of the total Virgo Gluster lumi¬ 
nosity. SMBBHs in low luminosity ellipticals contribute 
slightly more per galaxy to the ICL (3%) than the bright¬ 
est ellipticals, but most SMBBH ejections occur in the 
brightest magnitude bins, which are occupied by large 
ellipticals like M87. The contribution from spirals is 
negligible (0.01%). Of course, the fraction of light from 
galaxy interactions is more substantial; N-body simula¬ 
tions suggest 10 — 50% of the luminosity of the Virgo 
Cluster could come from intracluster stars stripped dur¬ 
ing galaxy interactions (Napolitano et al. 2003). How¬ 
ever, this study indicates that as much as 20% of the ICL 
could come from the more exotic SMBBH origin. We 
caution that since this estimate clearly hinges on M,, 
varying the M, — a relation can change Mej (and hence, 
the fraction of SMBBH-generated ICL) significantly. 

5. IMPLICATIONS 

Intracluster stars may be released via two very distinct 
mechanisms: tidal stripping of the outskirts of galaxies, 
and 3-body interactions with a SMBBH deep within a 
galaxy potential well. Stars ejected by a SMBBH are 
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kinematically separable from the ambient cluster kine¬ 
matics and also preserve an easily discernible shape, even 
in a complex and evolving cluster potential. 

A thin torus is the prompt ejection signal from a single 
SMBBH as it forms a hard binary. Hence, observing a 
torus of ICPNe can constrain both the recent existence 
of a SMBBH and its orbital plane without sub-arcsecond 
observations. For example, given a full loss cone, it takes 
~ 10^ years for a 10 : 1 SMBBH to evolve from a hard bi¬ 
nary to one that decays via gravitational radiation. Dur¬ 
ing that time, obeh shrinks, lowering the stellar interac¬ 
tion cross-section but strongly increasing the few ejected 
velocities. So, if the hypervelocity ICPNe observed are 
indeed caused by an interaction with a SMBBH, the last 
major galaxy merger was complete less than 10^ years 
ago. With a 10 : 1 galaxy merger timescale of 0(10®) 
years, there should be some fine structure (Schweizer & 
Seitzer 1992) within M87 associated with the event (and 
there may be: see Weil et al. 1997). 

The current consensus is that most ellipticals have as¬ 
sembled by merging (e.g., van Dokkum 2005). Therefore, 
in a cluster environment, debris may originate from sev¬ 
eral galaxies over many epochs, which makes the ejection 
signature more complex. To determine the aggregate be¬ 
havior of SMBBH-driven ICL we plan to investigate the 
long-term time evolution of the ejecta from SMBBHs as 
galaxies assemble within a cluster environment. 

The ICL fraction generated by SMBBHs (4 — 20%) 
is consistent with the high velocity ICPNe fraction in 
Virgo, 7%. Hence, the three ICPNe outliers could all 
be hypervelocity ejections. As a further reality check, 
given the predicted 0(10®) km/sec velocity kick, it takes 
O(IO^) years to travel 65 kpc from M87. Since the life¬ 
time of a 2Mq star is 0(10®) years, it is plausible that 
a main sequence star could plunge from the inner kpc, 
get ejected by the SMBBH, and evolve into a PNe within 
the allotted time. The tidal disruption radius of such a 
star is 0(10“®) pc, so it may pass within obbh intact. 

There is another way to determine the importance of 
SMBBH ejections. Since this mechanism works deep 
within the galaxy potential, the ejected stars should re¬ 
tain the metallicity of the elliptical galaxy center. Obser¬ 
vations indicate that bright ellipticals have a mean abun¬ 
dance in the range 0 < \_Z/H\ < 0.4 and a metallicity 
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gradient of AlogZ/Alogi? ~ —0.25 (Henry & Worthey 
1999). This implies that SMBBH ejection candidates 
will typically have a much higher metallicity than the 
tidally generated ICL, which usually originate in the out¬ 
skirts of a galaxy. It would be ideal if this metallicity di¬ 
vide were clear-cut, but some scatter in the abundances 
of both types of ICL is expected. Recall that, by far, 
most of the SMBBH ejecta originate from the center of 
the largest ellipticals, and should therefore exhibit high 
metallicities. This high metallicity sample is contami¬ 
nated, though, with a few low metallicity SMBBH ejec¬ 
tions that originate from faintest ellipticals (Bender et 
al. 1993). Strong galaxy interactions, too, can occa¬ 
sionally dredge up metal-rich material from the inner re¬ 
gions of galaxies (Hibbard & Mihos 1995; Mihos 2006 in 
prep), after which the cluster potential can easily strip 
it into the ICL. Finally, since the scatter in \Z/H]^ is 
large within a galaxy, metal poor stars exist even within 
a metal rich bulge (e.g., Jacoby & Ciardullo 1999). De¬ 
spite these complications, since the mean metallicities of 
each population differ, multi-night spectroscopic obser¬ 
vations with an 8-meter telescope (or better yet, short 
observations with a 30-meter class telescope) can divide 
the population statistically. 

The implications for LISA are profound. If a star 
passes within 100 Schwarzschild radii of either SMBH, 
it can generate a burst of gravitational radiation with 
a signal-to-noise ratio of 0(10®) in the LISA waveband, 
high enough to be observed in Virgo at a rate of 10 per 
year (Rubbo, Holley-Bockelmann, & Finn, in prep.). In 
the future, we may use LISA observations in concert with 
hypervelocity PNe to confirm and further constrain the 
masses, mass ratios, separation, and even the SMBBH 
orbital plane within galactic nucleii. 
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